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Abstract 

Negative ions (anions) were identified in the coma of comet IP/Halley from in-silu measurements 
performed by the Giotto spacecraft in 1986. These anions were detected with masses in the range 7-110 
amu, but with insufficient mass resolution to permit unambiguous identification. We present details of 
a new chemical-hydrodynamic model for the coma of comet Halley that includes - for the first time - 
atomic and molecular anions, in addition to a comprehensive hydrocarbon chemistiy. Anion number 
densities are calculated as a function of radius in the coma, and compared with the Giotto results. 
Important anion production mechanisms arc found to include radiative electron attachment, polar 
photodissociation, dissociative electron attachment, and proton transfer. The polyyne anions C 4 H' and 
Cf,!-!' are found to be likely candidates to explain the Giotto anion mass spectrum in the range 49-73 
amu. The CN' anion probably makes a significant contribution to the mass spectrum at 26 amu. Larger 
carbon-chain anions such as C^H' can explain the peak near 100 amu provided there is a source of large 
carbon-chain-bearing molecules from the cometaiy nucleus. 


1. Introduction 

Gas-phase negative ions have been detected in surprisingly large abundances during spacecraft fly-bys 
of Solar System bodies. During the March 1 986 encounter of the Giotto spacecraft with the coma of 
comet IP/Halley, the RPA-1 Electron Electrostatic Analyzer (EEA) detected energy spectra consistent 
with the presence of cold anions in the coma with masses 7-110 amu and number densities up to 1 cm'^ 
(Chaizy et al. 1991). Atomic and molecular anions O', OH', CN’, C^IT and larger CHON' particles 
were proposed by Chaizy et al. ( 1 99 1 ) as plausible candidates to explain the observations, but due to 
the limited mass-resolution of the Giotto anion mass spectra, unambiguous identification of the 
detected anions was not possible. The Cassini spacecraft performed similar measurements using the 
CAPS Electron Spectrometer (ELS) within the atmospheres of Titan (Coates et al. 2007) and Bnceladus 
(Coates et al. 2010), and detected negative ions with masses per unit charge ranging from ~ 10 to 
1 0,000 amu/q. 

Molecular anions were first detected in the interstellar medium using microwave spectroscopy 
(McCarthy et al. 2006). They have since been found to be abundant in a range of environments outside 
the Solar System, and are an important part of the molecular inventory of the Galaxy (Cordincr & 
Chariiley 2012). Carbon-chain anions have been found to be abundant in quiescent molecular clouds, 
prestcilar cores and protostars (Brtinken et al. 2008, Gupta et al. 2009, Sakai et al. 2007, 2010, Cordiner 
et al. 2011), as well as the carbon-rich AGB star IRC-t-10216 (Remijan et al. 2007). Chemical models 
have been able to reproduce the observed anion abundances in interstellar and circumstellar 
environments, and have shown that molecular anions have important effects on the chemistry and 
ionization balance of the gas (Cordiner & Millar 2009, Walsh et al. 2009, Cordiner & Charnlcy 2012, 
McElroy et al., Forthcoming). Whereas chemical kinetic models have been successful in explaining the 
abundances of small molecular anions in Titan’s upper atmosphere (Vuitton et al. 2009), the origin and 
importance of anions in the coma of comet IP/Halley, as well as for comets in general (Wekhoff 1981), 



remains to be understood. 


In this article, we utilize a combined chemical/liydrodynamic model for the coma of comet Halley to 
explore various anion production mechanisms and compute the abundances of atomic and molecular 
anions as a function of radius in the coma. 


2 , Chemical model 

Our model is based on the five-fluid coma model of Rodgers & Charnley (2002), which models the 
dynamics of a neutral fluid, a positively-charged fluid, an electron fluid and fluids of fast atomic and 
molecular hydrogen. These fluids emanate from the nucleus in a spherically-symmetric outflow (see 
also Rodgers ct al. 2004), and the abundances of various chemical species arc calculated as a function 
of radius. The chemical network has been fully updated in the present work and incorporates reactions 
between 279 chemical species (composed of H, C, N and O atoms, and electrons), The list of species 
consists of 154 cations, 23 anions (listed in Table 1) and 101 neutrals. Hydrocarbon chemistry is 
modeled for species containing up to eight C-atoms. Carbon-chain-bcaring anions arc included because 
of the importance of C 2 and C 3 in cometary comae (Weiler 2012), and because their gas-phase 
chemistiy is better understood than other classes of anions such as PAHs. The abundances of the 
modeled species are linked by a total of 3823 chemical reactions, 3367 of which were taken from the 
latest version of the UMIST database for astrochemistry (RATE12; McElroy et al. Forthcoming), 
including 685 reactions involving anions. 


Table 1. Anions included in the model 


ir, O', OH', 02 *, 

C, C2', C3', CT, CT, CT, c/, CT, 

CH‘, C2H', C^ir, C4H', CsH', CftH', C7H', CgH', 
CN', C3N', C^N' 


Solar photo-reaction rates at 1 AU have been taken from Huebner et al (1992) and Crovisier (1994) 
where available. New photodissociation rates were calculated for hydrocarbons based on the cross 
sections of van Hcmcrt & van Dishoeck (2008), integrated over the quiet sun solar irradiance spectrum 
(Woods et al. 2009). For the remaining neutral species, known photo-rates for structurally similar 
species were used, with product channels taken from the RATE 12 database. Anion photodetachment 
rates were calculated using the square-root threshold law (Millar et al. 2007). We used measured cross 
sections above threshold for C 2 H', C^H' and C 6 H' (from Best et al. 201 1), and for OH" (Trippel et al. 
2006). Above-threshold photodctachmcnt cross-sections were assumed to be 10"'^ cm^ for all other 
anions. 

Mutual neutralization reactions have been included between all anions in the model and the most 
abundant cations, including H 20 ^, H 3 O", NHT, C 2 H 2 ", C 2 HT, HCNH^', CH.iOH.^ and CH 3 OCH/, with 
rate coefficients of 7.5 x 10’^(7’,-/300)'^ ‘'’cnr^s'*. In this formula 71- is the effective collision 
temperature, defined by 



j _ 

' 

where 7} and Tt, are the respective temperatures of the (cat)ionic and anionic reactants, and m, and m., 
arc their respective masses. 

A further 78 reactions involving anions have been added, including proton transfer, dissociative 
electron attachment and polar photodissociation. A summary of the main chemical processes involving 
anions in the coma model is given in Table 2. The dominant anion formation processes arc described in 
more detail in Section 3. 


Table 2. Anion reaction types included in the coma model 


Abbreviation 

Reaction Type 

Example 

REA 

Radiative electron attachment 

A + C A" + hv 

PPD 

Polar photodissociation 

AB + hv A' + B ■ 

DEA 

Dissociative electron attachment 

AB + C ^ A" + B 

PT 

Proton transfer 

AH + B*^ A' + BH 

MN 

Mutual neutralization 

A + B A + B 

AED 

Associative electron detachment 

A' 4 - B AB + e 

AN 

Anion-neutral 

A' + B ^ C + D 

CT 

Charge transfer 

A' + B A + B' 

PD 

Photodetachment 

A' + hv A + e' 


The contribution of anion chemistry to the energetics of the gas is negligible in the majority of models 
considered here due to the relatively small anion abundances. The kinetic temperature of the anions is 
assumed to be identical to the temperature of the neutral gas from which they arc formed. 

Our model for Halley assumes a spherical nucleus with radius 3.36 km, heliocentric distance 0.89 AU, 
temperature 192 K and H 2 O production rate of QflhO) 6.9 x 10"^ s‘‘, appropriate to the time of the 
Gioffo encounter on 1 3 March 1986 (Schmidt ct al. 1 988). Production rates of parent species are given 
in Table 3. These have been taken from the compilation of Haider & Bhardwaj (2005), with the 
addition of C«H 2 , that has been included as representative of a generic long-chain hydrocarbon 
molecule (see Section 5). 


Table 3, Production rates for parent species relative to II 2 O 


Species 

Production 

Rate 

Species 

Production 

Rate 

H2O 

1.0 

C2H4 

3.00 X 10-’ 

CH4 

6.25 X 10-’ 

CO 

1.00 X 10‘‘ 

NIC 

1.88 X 10’^ 

C2H, 

4.00 X lO'-'' 

C2H2 

2.50 X 10'- 

CH3OH 

2.13 X 10'“ 

HCN 

1.25 X I0‘" 

CO2 

3.75 X 10'“ 

N2 

1.25 X 10" 

CgH. 

1.00 X 10" 





3. Anion production processes 


3.1 Radiative electron attachment 

Two main processes are responsible for initiation of the anion chemistry in the inner coma (at radial 
distances leSvS than ~ 1000 km from the nucleus); radiative electron attachment and polar 
photodissociation. As a result of their large electron affinities, combined with a high vibrational 
density of states, carbon-chain-bcaring molecules (including polyynes, C„H, for n > 3), possess the 
ability to attach low-energy free electrons and undergo rapid radiative relaxation to form stable anions 
(Herbst & Osamura 2008). As sliown in Fig. Ic, the Cf^H' number density quickly reaches > 1 cm'"* as a 
result of radiative electron attachment to CisH. The C(sH molecule is formed in the coma from electron 
recombination of Cc,H 4 ^ - a moderately large hydrocarbon ion that arises in the model as a product of 
acetylene photochemistry. The neutral-neutral reaction 


C2 + C2II2 C4H + H 


results in the formation of a large amount of C4H (reaching a production rate of ~ 10'“^ (2(^^20)). 
However, due to its smaller size, C4H is theorized to possess a smaller radiative electron attachment 
rate than CrJl (Herbst & Osamura 2008), causing C 4 H' to be less abundant than C6H' in the inner eonia. 
The larger CglT anion is abundant in our model as a consequence of the inclusion of Cgl H as a parent. 
In the innermost regions of the coma (at radial distances <10 km), where few Solar UV photons 
penetrate, CgH is produced as a result of charge transfer from H" to C 8 H 2 , which forms C 8 H 2 ~ ions that 
undergo subsequent electron recombination to CgH + H. At larger radii, photodissociation of Csl-h is 
the dominant source of CgH. 

The majority of the C 3 N' in the model comes from radiative electron attachment to C 3 N, which is 
produced predominantly as a result of the neutral-neutral reaction 


C 2 + HCN ^ C 3 N + H. 


3.2 Polar photodissodation 

At distances greater than ~ 10 km from the nucleus, the penetration of solar UV radiation becomes 
sufficiently strong that polar photodissociation (PPD) of neutral molecules becomes an important 
source of negative ions. Also known as threshold ion pair production, PPD is the process by which 
absorption of radiation by a molecule AB results in the production of a pair of oppositely-charged 
product species A' and B' (sec Table 2). Polar photodissociation of H 2 O, HCN, CO and CO 2 was 
studied in the laboratory by Hunniford ct al. (2007), Berkowitzet al. (1969), Dadouch (1991) and 
Mitsukc (1990), respectively. Our model incorporates rates for PPD calculated in the solar radiation 
field using the cross sections given by Vuitton et al. (2009), integrated over the applicable energy range 
for each process (taken from the original laboratory studies). 
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Fig. 1. Modeled negative ion (anion) number d< 
distance from the nucleiKs. 






Polar photodissociation of H 2 O by photons in the energy range 36-100 cV produces both O' and H' 
(Hunniford et al. 2007). This is the dominant production mechanism for atomic anions in our 
coma model. The O' number density as a function of radius is shown in Fig. la, which reaches a peak 
value of ~ 1 cm'"' at a radius of - 80 km. Polar photodissociation of CO is the primary production 
mechanism for C throughout the coma. 


3,3 Proton transfer 

Polar photodissociation of HCN is an important source of CN', but the dominant production mechanism 
for this anion in our model is via loss of protons from HCN through proton-transfer reactions. The 
transfer of protons from neutral species to anions can be highly cxoergic depending on the acidity of 
the neutral in question, and such reactions are almost always extremely fast in the gas phase (Vuillon et 
al. 2009). For the anions in our model, we have included all relevant exocrgic proton transfer reactions, 
with rate coefficients taken from the RATE 12 database and from Vuitton el al. (2009) where available. 
Reaction rates have been extrapolated (o the additional anions in our model where necessary. 

As a result of its high proton affinity, H' quickly reacts with l EO to produce OH' + H 2 . Consequently, 
OH' reaches a large number density in the coma, similar to that of O'. Abundant C 2 H' is produced by 
proton transfer from C 2 H 2 to O' and OH'. Proton transfer from HCN to C4H' and C(,H' is the dominant 
source of CN' for radii less than ~ 2,000 km. 


3.4 Dissociative electron attachment 

For radii greater than ~ 2,000 km, dissociative electron attachment (DEA) to HCN is the dominant 
source of CN'. The onset of this formation process manifests as the emergence of a secondary peak in 
the CN' number density al around 10,000 km (sec Fig. lb), dominating over the proton transfer and 
polar photodissociation formation channels at these radii. Dissociative electron attachment rates have 
been calculated for the anions in our model by convolving the (Maxwellian) electron energy 
distribution with experimentally-derived cross sections for anion formation as a function of energy. 
References to much of the original laboratoiy data are summarized by Vuitton et al. (2009), who also 
provide the peak cross sections and energies. Additional DEA cross section data for C 2 II 2 , CO 2 and O 2 
have been taken from May et al. (2008) and Rapp et al. (1965). Gaussian profiles have been assumed 
for the cross sections as a function of energy, which is a good approximation in most cases. A Gaussian 
FWHM of 1 .0 cV has been used in the case of HC.-?N (and HC 5 N), for which no detailed cross section 
profile data have been published. 

Dissociative electron attachment only becomes significant as a source of negative ions in the outer 
coma, where the electron energies become sufficient to overcome the energy barrier of this 
endothermic process. For example, DEA of HCN occurs at an electron energy of 2.5 cV, so the peak 
rate for CN' formation from HCN occurs for electron temperatures Te ~ 2.9 x 10“^ K, but begins to 
become significant at > 1000 K, when sufficient electrons possess the required 2.5 eV. Energetic 
electrons arc produced by photo-ionization of water molecules. However, their temperatures only 
obtain large values in the outer coma where the lower density results in a reduction in the amount of 
energy lost in inelastic collisions with H 2 O (Rodgers & Charnlcy 2002). 

Secondary anion abundance peaks in the outer coma are also observed in our model for the following 
anions (see Fig, 1); OIF (produced from DEA of H 2 O), C 2 H', C 2 ~ (from DEA of C 2 H 2 ), and C 3 N' (from 



DEA of HC3N). Secondary DEA peaks arc not seen for the larger polyync anions (C„ir, for n > 3), 
because their large radiative electron attachment rates dominate the other anion formation processes. 


5. Discussion 

Chaizy et al. (1991) detected anions in the .coma of comet Halley with masses in the range 7-110 amu, 
and despite the limited resolution of the mass spectra, identified three main peaks at (1) 7-19 amu, (2) 
22-65 amu and (3) 85-1 10 amu. Tliey referred to these as the ‘17 amu group’, the ‘30 amu group’ and 
the ‘100 amu group’, respectively. The most likely candidate species to explain Peak 1 are shown in 
Fig. la: O' and OH', C and CH'. In our model, O' is the most abundant of these anions throughout 
much of the coma - OH' is slightly more abundant in the inner regions (at a radial distance r < 200 km 
from the nucleus) and outer regions (r > 10,000 km). Chaizy ct al. (1991) derived anion number 
densities for Peak 1 starting at more than 1 cm"'^ at ;• = 3 x 10'^ km and falling to ~ 10'"^ cm'^ at r - 3 x 
10"^ km, (as shown in their Fig. 4). Outside of the comet’s ionopause (at r = 4,500), our modeled Peak 
1 anion abundances are about two orders of magnitude less than observed. 

The mass range covered by Peak 2 contains the relatively small di/tri-alomic anions C 2 ', CN‘ and C 2 H', 
as well as some of the medium-sized carbon-chain-bearing anions: C4H' and C 3 N'. Over the range of 
the Giotto data (from r = 3 x 10’^ km to 3 x 10"^ km), the observed Peak 2 anion densities fall from ~ 10' 
■ to 1 O'" cm'\ As can be seen in Fig. I , our modeled number densities for CN' and C 4 H‘ match this 
trend quite well (within an order of magnitude) al all radii. Fig. 4 of Chaizy ct al. (1991) shows a rise in 
the observed number density of these ‘30 amu group’ anions between 1 x 10"^ km and 2x10'* km, 
which is consistent with the location of the rise in anion abundances (including that for CN'), which 
occurs as a result of DEA in the outer coma. 

Peak 3 contains the larger carbon-chain-bcaring anions CfiH' and CglT. Our modeled number densities 
for these species as a function of radius are again consistent (within an order of magnitude) with the 
number densities of anions in this mass range detected by Giotto. 

The CglE molecule is included as a parent in the coma model, and is intended to be representative of 
the class of large carbon-chain-bcaring species that could plausibly be present as a component of 
conietary nuclei (Mumma & Charnlcy 201 1). Photodissociation of C^IH is assumed to result in H-atoni 
loss to produce the polyyne CgH, although alternative product channels are possible, including loss of 
various numbers of carbon atoms. Long carbon chains and polyynes (including CgH) undergo radiative 
electron attachment at similarly rapid rates (Herbst & Osamura 2008). Thus, CgH' in our model can be 
considered as representative of the class of anions that arise in the coma as a consequence of the release 
of dusty, carbon-chain-rich material from the nucleus. The appearance of a strong anion mass peak in 
the Giotto observations at ~ 100 amu is thus readily explainable. As a final point, we note that if a 
sufficient quantity of long carbon chains is released from the nucleus (with a production rate greater 
than 1 0 '*^ 0 (H 2 O)), rapid radiative electron attachment ensues and molecular anions can be formed in 
such large abundances that they become the dominant carriers of negative charge in the inner coma (at 
distances < 10'^ km from the nucleus). Evidence for the presence of abundant carbon chains in Halley 
was provided by Geiss et al. (1999), who derived a C4H production rate of 0.023 ^(ILO). In such a 
situation, the total abundance of anions is sufficiently large that their thermodynamic effects on the 
coma begin to become important. A proper treatment of these effects will require the inclusion of a 
separate anion fluid in future hydrodynamical coma models. The properties of such an anion fluid will 
be examined in detail by Cordiner & Charnley (Forthcoming). 



The scatter on the anion number densities as a function of radius presented by Chaizy ct al. (1991) is 
considerable, and often larger than the associated error bars. Uncertainty is inherent in their results due 
to the fact that the Giotto EESA was calibrated for use in electron spectroscopy rather than (heavy) ion 
spectroscopy. In fact, the efficiency of the EESA microchannci plate (MCP) detector for measuring 
negative ion fluxes is highly uncertain due to a lack of faboratoiy calibration for such purposes. Thus, 
the results of Chaizy et al. ( 1 99 1 ) can only be considered approximate such that a direct quantitative 
comparison of our results may not be meaningful. Nevertheless, the qualitative agreement between 
observations and our model predictions is good, especially considering some of the uncertainties in the 
model. Based on laboratory measurements of polyyne anions and OH‘, the photodetachment cross 
sections employed for many of the other anions (including O', OH', CN' and CjN'), are likely to be 
accurate to within a factor of a few al best. Because photodetachment is by far the dominant anion 
destruction mechanism, errors on these cross sections propagate in an approximately linear fashion to 
the modeled abundances of these species. For those species whose production is dominated by polar 
photodissociation, a potentially more significant source of error is the uncertainty in the PPD cross 
sections, which have been approximated and may only be accurate to within 1-2 orders of magnitude. 
Thus, we highlight a need for dedicated laboratory measurements of these values. Due to a similar lack 
of laboratory measurements for electron attachment to neutrals at low energies and densities, the 
radiative electron attachment (REA) rates for all anions in the model are based on theoretical 
calculations. Previous endeavors to model observed anion abundances in dark molecular clouds and 
circumstcllar envelopes have shown that the theoretical REA rates are reasonably accurate for CgH' and 
CsH', but less so for C 4 IT, the rate for which seems to have been over-estimated. 

Other possible sources of uncertainty in our model include a lack of measured rate coefficients for 
several important proton transfer reactions, including the reactions of O' and OH' with C 2 H 2 , and the 
polyyne anions with HCN. In addition, our list of PPD and DEA reactions is incomplete due to a lack 
of laboratory studies of these processes for some abundant coma molecules such as NH 3 , H 2 CO and 
CH 3 OH. Possible anionic products from PPD and DEA reactions involving these species may include 
NH 2 ', CH' and CII 30 ', respectively, whieh are not included in our model. 

The fact that our model produces about two orders of magnitude less Peak 1 anions (O', OH') than 
observed may be indicative of missing formation processes for these species. For example, Chaizy ct 
al. (1991) hypothesized that anions could be produced in high-energy collisions between neutral 
species in the coma. The process of anion production through collisions between neutral species at 
energies E up to 20 eV was reviewed by Wexler (1973). Cross sections for anion production in higher- 
energy impacts (£■- 100 - 1,000 eV) are typically very small (McDaniel 1964), but a value of ~ 10''^ 
cm” at 100-250 eV was measured by Gealy & van Zyl (1987) for anion production in collisions 
between H and H 2 . This process could be an important source of anions in the coma given a sufficiently 
strong source of energetic atoms or molecules. Eviaiar et al. (1989) deduced that such a flux of neutral 
species may be present in the coma of comet Halley in order to explain the fluxes of fast ions detected 
by Giotto. Given the dominance of water in the coma, collisional dissociation of H 2 O may thus be 
considered a plausible source of H', O' and OH'. However, collisional dissociation of H 2 O - and other 
abundant coma molecules - has yet to be studied in the laboratory, so the existence of these 
hypothetical anion production channels (not to mention their cross sections), is presently unknown. 



5. Conclusion 


A reasonably good agreement has been achieved between observed and modeled anion abundances in 
the coma of comet 1 P/Halley, despite uncertainties in various cross sections and rates for anion 
production and destruction. The dominant anion production mechanisms arc found to be polar 
photodissociation of water and radiative electron attachment to carbon chains in the inner coma, 
followed by proton transfer from C 2 H 2 and I ICN to produce C 2 H’ and CN', respectively, In the outer 
regions of the coma where electron temperatures reach ~ 1 0^-10’ K, dissociative electron attachment 
becomes a dominant process, 

We find particularly good agreement for CN' and C 4 H' in the ‘30 amu group’, and for C^jH' and C^H' in 
the ‘ 1 00 amu group’. We thus confirm the hypotlicsis that CN' is a likely carrier of the 20-30 amu 
anion mass peak observed by Gioito in the 1 986 encounter with comet Halley. The polyyne anions 
C 4 H' and Cc,H' arc likely constituents of the mass spectrum in the range 49-73 amu. Larger carbon- 
chain-bcaring anions such as C^' and C^H' can explain the mass peak near 100 amu, provided a source 
of large carbon-chain-bearing molecules is present in the cometary nucleus. 

If our interpretation of anion chemistry in the Giotto data is coirect, then comets may contain 
significant abundances of long carbon-chain molecules. Measurements with the Ion and Electron 
Sensor on the Rosetta spacecraft during its scheduled 2014 encounter with comet 67P/Churyumov- 
Gerasinienko will provide additional insight into this issue. A more complete understanding of 
cometary anions will be facilitated by future laboratory measurements of absolute cross-sections for 
anion photodetachment and polar photodissociation of H 2 O and other abundant coma molecules. 
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